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Abstract 

Recently we have described a mechanical system which exhibits 
spontaneous breaking of symmetry and related topological kinks 
called compactons. The corresponding field potential is not differen- 
tiable at its global minima. Therefore, standard derivation of disper- 
sion relation u(k) for small perturbations around the ground state can 
not be applied. In the present paper we obtain the dispersion relation. 
It turns out that evolution equation remains nonlinear even for arbi- 
trarily small perturbations. The shape of the resulting running wave is 
piecewise combined from ±cosh functions. We also analyze dynamics 
of the symmetry breaking transition. It turns out that the number 
of produced compacton-anticompacton pairs strongly depends on the 
form of initial perturbation of the unstable former ground state. 
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1 Introduction 



During the recent decades it has become clear that spontaneous symmetry 
breaking and related with it topological defects play important role in physics 
of superfluid Helium, superconductors, liquid crystals, hadronic matter, cos- 
mic objects like neutron stars or hypothetical cosmic strings, and many other 
systems. Reviews and references can be found in, e.g., [HE]. Static prop- 
erties of the topological defects are understood relatively well. Dynamics of 
them is much harder to unravel. Interactions and production of the defects 
in general involve many modes of underlying fields (or order parameters), 
and non-trivial time-dependent solutions of nonlinear field equations have to 
be investigated. 

In the present paper we discuss time-dependent solutions in the effective 
model introduced in The model is distinguished by several attractive 
features. First, the pertinent field equations are very simple. Second, the 
model has clear relation with a classical system of macroscopic pendulums, 
which is easy to visualize, and even to build. This is in contrast with many 
other effective models with topological defects where their relation with the 
pertinent physical system is not transparent - this is especially true in the 
case of classical Ginzburg - Landau type effective models for quantum fluids. 
Finally, the topological defects in our model have only contact interactions - 
they do not interact until they touch each other. The reason is that solutions 
describing the defects exactly coincide with vacuum solutions outside certain 
finite region in the space, hence the common exponential tails are absent 
here. Such defects are called compactons 1 . 

The Lagrangian obtained in [2j has the form 

L = l(d T <j ) ) 2 - 1 -(d^r-v(<j ) ), a) 

where £ is a dimensionless coordinate along a straight-line, = 0(£, r) is the 
scalar field which in general can depend on £ and on a dimensionless time r. 
The field potential V(<p) is given by the following formula 

J cos0-l for |0| < O 
VYV [ oo for |0| > 0o, w 

where 0o is a constant, n > 0o > 0. The potential has two degenerate min- 
ima at = ±0o, see Fig. 1. 



1 This name has been proposed in connection with solitonic solutions of modified Ko- 
rteveg - deVries equations, see, e.g., pHO^E). 
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V(<|>) 




Fig.l. The potential V(4>). 



The variable £ is restricted to the interval [— 6v,6v]. By assumption, the 
field obeys the following boundary conditions 



Because of the infinite potential barrier, modulus of the field cannot exceed 

00, 



Scalar fields with values restricted to certain finite interval appear also in a 
completely different context of brane models, [Zj. 

As discussed in detail in [3j, the effective model defined by formulas (1-4) 
describes a system of 2N + 1 pendulums attached to a rectilinear wire at the 
constant distance a between them. Each pendulum has the arm of length R, 
and mass m on the free end which can move only in a plane perpendicular to 
the wire. is the angle between the arm and the vertical upward direction. 
The pendulums are firmly attached to the wire, so that is also the angle by 
which the wire is twisted. It is assumed that the wire is elastic with respect 
to the twists, k is the corresponding torsional elasticity coefficient. The angle 
is restricted to the interval [— 0o, 0o] by two stiff rods placed symmetrically 
on both sides of the wire, parallel to it and at the appropriate distance. The 
boundary conditions (3) mean that the two outermost pendulums are kept at 
rest in the upward position by an external force. The dimensionless variables 
£, r are related to the physical position coordinate x (along the wire) and 
time t by the formulas 



0(±6v,r) = O. 



(3) 



- 00 < < 00- 



(4) 
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The two outermost pendulums are located at x — ±Na, hence 



(6) 



The system of pendulums as described above, but without the two rods, 
is well-known in connection with sine-Gordon equation and its solitonic so- 
lutions [HJ. In that case the space of values of r) is a circle S 1 , and the 
solitons appear because of nontrivial topology of S 1 . In our case, instead 
of S" 1 we have the topologically trivial interval [— 4> ,(p }. The presence of 
compactons is due to spontaneous breaking of the Z 2 symmetry 



It takes place when £jv becomes large enough, see Section 3. Yet another sys- 
tem of coupled pendulums is described in jHj. It differs from the sine-Gordon 
or our system by the presence of a term which represents nonlinear elasticity 
of the wire, and from our system also by absence of the rods. Thanks to the 
nonlinear elasticity the solitons are compactons. They essentially differ from 
the ones discussed by us. 

In |H1 we have discussed in detail the relation between the system of 
pendulums and the field-theoretic model with Lagrangian (1), as well as basic 
properties of compactons: their shape and energy. In the present paper we 
investigate waves propagating in the ground state with broken Z 2 symmetry, 
and we make a preliminary study of production of compactons and anti- 
compactons during the Z 2 symmetry breaking transition. Such transition 
can be triggered by, for example, decreasing the elasticity coefficient k or 
by increasing the length R of the arm. Both possibilities seem technically 
feasible. The first topic is discussed in Section 2. It is interesting because 
the standard procedure, based on linearization of field equations, does not 
work here - potential (2) is not smooth at the absolute minima. It turns out 
that even in the limit of small amplitudes the pertinent evolution equation 
remains nonlinear. As for the second topic, it is related to the broad subject 
of dynamics of rapid phase transitions and production of topological defects 
(KB El- We show in Section 3 that thanks to extraordinary simplicity of our 
model one can calculate time evolution of the system during the symmetry 
breaking transition. In particular, one can see how many topological defects 
are produced. 

2 Small perturbations of the ground state 

In this section we consider an infinite chain of pendulums, — * an d 
we abandon the boundary conditions (3). The main reason for this step is 



0(£,t). 



(7) 
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that we would like to study effects which are due to the infinite potential 
barrier at = ±0 O . At the ends of the finite chain the barrier is not felt at 
all because of condition (3). 

The evolution equation which corresponds to Lagrangian (1) has the form 

^0-d|0-sin0 = O (8) 

if 

101 < 00- 

When reaches ±0o the pendulums elastically bounce from the rods, 

<9 T - -d T <t>, (9) 

or just rest on them if the velocity vanishes. Let us emphasize that Eq.(8) 
does not hold when = O . 

In the present paper we assume that the maximal angle 0o is small and 
therefore 

sin0 = 0. (10) 

This can be achieved by putting the rods close enough to the pendulums 
when they are in the upward vertical position. Then, Eq.(8) is linearized 

<^0 - df0 - = 0. (11) 

For the infinite chain of pendulums the two degenerate ground states are 
given by 

= ±00- (12) 

Evolution of small perturbations of the ground state 0o is governed by Eq.(ll) 
and the reflection condition (9). Such perturbations have the form of waves 
propagating along the chain: the pendulums rise a little bit above the rod 
and fall back. We assume that amplitudes of the waves, equal to 0o — 0(C) r )i 
are not too large, 

O - 0(£, r) < O , (13) 

that is that none of the pendulums reaches the upward vertical position. 
Condition (13) is of course equivalent to 0(£, r) > 0. In order to investigate 
the waves we use the following trick. Instead of the original field 0(£, r) we 
consider another one, namely 0(£, r) related to 0(£, r) by the formula 

™> T) \ 20 o -0(£,r) if 0(£,r)>0 o . [M) 
Condition (13) means that 

< < 20 o . (15) 
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Equation (11) is equivalent to the following evolution equation for 



(a? - = 



if < 00, 

if = 00 

0-200 if > 00- 



(16) 



The r.h.s. of this equation for ^ O corresponds to the potential K(0) 
given by the following formula 



Z(0) = 



for < < 



|0 2 iur u <^ (p \ (/^o, 

■|(0-20 o ) 2 for 0o<0<20 o . 



(17) 



H(0) is symmetric with respect to the point = 0o at which it reaches its 
absolute minimum, V_{4>o) = — 0o/2, see Fig. 2. Equation (16) gives contin- 
uous evolution of and of its first derivatives <9 T 0, <9g0. The advantage of 
using instead of is that because of relation (14) such continuous evolution 
of <9 T automatically gives the reflection (9) of <9 T at = O . 




Fig.2. The potential V(f). 
Because we are interested in oscillations of around O , we substitute 

= 0o + o e(£,r). 
Evolution equation for e follows from Eq.(16), 

(d 2 T - %)e = e - sign(e) (18) 
(we assume that sign(O) = 0). Condition (15) is equivalent to 

|e| < I. (19) 
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Solutions of Eq.(18) describing running periodic waves can be obtained with 
the help of the Ansatz 

e(£,T) = e(C), C = 4^=- ( 20 ) 
V v 2 — 1 

Here v is the phase velocity of the wave. It is assumed that |i>| > 1. Now 
equation (18) is reduced to ordinary differential equation 

e" = e - sign(e), (21) 

where ' denotes d/d(. Integration of Eq.(21) gives the relation 

6' 2 = (l-| e |) 2 -(l-60) 2 , 

where eo > is the maximal value of e. Hence, 



±(? = y/e -\e\y/2-e -\e\. (22) 

the sign + or — should be chosen in accordance with the sign of e' because 
the r.h.s. of (22) is nonnegative. We consider Eq.(22) in the regions e > 
and e < separately. Simple calculations show that if ±e > then e = e±, 
respectively, where 

e+(C) = l-(l-eo)cosh(C-Co), (23) 

e_(C) = -l + (l-eo)cosh(C-Ci)- (24) 
Here Co, Ci are constants. These two solutions should match each other at 
C — Ci/2 such that e + = = e_. Thus, 

1 = (1 - eo) cosh(Ci/ 2 - Co), 

and 

1 = (1 - eo) cosh(Ci/ 2 - CO- 
These conditions are satisfied if 

Ci/2 = Co + ^, Ci = Co + ^, (25) 

where 

A = 4 arcosh- . (26) 

1 — eo 

It turns out that at such C1/2 also the first derivatives of e + ,e_ are equal. 
e+(C) i s non-negative in the interval ( G [Co — A/4, Co + A/4], and e_(C) is 
non-positive on the interval C G [Ci — A/4, ^ + A/4]. Taken together they 
give e in the interval [Co — Co + |A], and e for all other values of C can be 
obtained by periodic repeating e + and e_, see Fig. 3. Thus, we have obtained 
the running periodic wave of the cosh type. 
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Fig. 3. The periodic wave of the cosh type. 



It is clear that A is the wave length for e(£). The periodicity in ( implies 
periodicity in £ with the wave length 




and in r with the period 



The corresponding wave number k = 2n/\ and the frequency u> = 2n/T 
obey the relation 

4vr 2 

- 2 - k 2 = %. (27) 
The group velocity of this wave 

duo k 1 
9 dk uo v 

is smaller that 1. Because the r.h.s. of formula (27) does not vanish, the 
waves are massive. For fixed wave number k the group velocity increases 
with A. According to formula (26) A increases with e . Therefore, the waves 
with larger amplitude are faster. This dependence on the amplitude appears 
because Eq.(21) is nonlinear. Due to the discontinuous character of the sign 
function it can not be linearized even for arbitrarily small e - in our system 
even arbitrarily small oscillations around the ground state are nonlinear. 



3 The Z<i symmetry breaking transition 

In this Section we consider the finite chain of pendulums, £at < oo. obeys 
the boundary conditions (3). We again assume that 0o is small enough to 
justify the approximation (10). In order to simplify the model further, we 
now assume that the dynamics of pendulums is dissipative, that is that the 
pertinent evolution equation has the form 

d T (f) - <9f -0 = 0, (28) 
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<t>v = < 



instead of (11) (if |0| < <f> .) Physically this can be achieved by putting the 
chain of pendulums in a viscous liquid. In the present case the second of 
formulas (6), relating the r variable with the physical time t, is replaced by 

mg 

where 7 is a constant coefficient characterizing the friction. The main ad- 
vantage of such over-damped system is that by quickly loosing its energy it 
rapidly evolves towards a stable state. 

If the chain of pendulums is long enough, that is if £jv > 7r/2 as shown 
below, there are two degenerate ground states given by 

= ±4>v, 

where 

O sin(£ + 6v) for £ G [-£ N , -£ N + §], 
0o for £g Kiv + f,&v-f], (29) 

_ 0o sin^Tv - for £ G - § , £ N ]. 

Such (f> v obeys the boundary condition (3). The sinus functions are solutions 
of Eq.(28) - they describe gradually tilted pendulums. For the intermediate 
values of £, that is when — £/v + n/2 < £ < £jv — 7r/2, the pendulums just rest 
on one rod. For the chains shorter than n formula for 0„ given above can 
not be accepted. Below we show that in this case the stable ground state is 
given by 

S (O = O. (30) 

This state is symmetric with respect to the Z 2 transformation (5). 

The static compacton is obtained from the ground state — <j) v by flipping 
pendulums from one rod to the other one at the intermediate values of £. 
This flip is given by the static solution of Eq.(28), 

= o sin(£-£o) (31) 

in the interval £ G (£0 _ 7r /2,£o + tt/2). The compacton is located at £0 
and it has the width n. Notice that the sinus functions present in formula 
(29) can be regarded as half-compacton and half-anticompacton sitting at 
the ends of the wire. If the compacton is present, at the ends we have 
two half-anticompactons. The compacton does not interact with these half- 
anticompactons if it is located far enough from the ends, namely if — £tv + 7t < 

£0 < 6v - 7T. 

Solutions of Eq.(28) can be written in the form of decomposition into 
eigenmodes 

00 Tin 
0(£^) = ]T/Jr) sin -^(£ + 6v), (32) 

n=l Z ^ N 
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where 



fn{r) = /„(()) exp(A n r), 



(33) 




(34) 



Let us recall that 0(£, r) evolves according to equation (28) until a pendulum 
touches one of the rods. Therefore, solution (32) is physically relevant until 
|0| = 4>o at certain £ and r. At later times the description of the evolution of 
the system has to take into account the reflection condition (9). Obviously, 
formulas (32-34) correctly describe evolution of small deviations from the 
Z 2 symmetric state S = 0. In the case when this state is unstable, these 
formulas are relevant if r is in certain finite interval [0,r ]. If that state is 
stable, A n < 0, the formulas describe exponential disappearance of small 
perturbations when r — > oo. Of course, in any case we have to take the 
coefficients /„(0) small enough. 

Formulas (33-34) show that the modes with n > 2£n/it exponentially 
disappear, while the modes with n < 2^n/tt exponentially grow. Because 
the smallest value of n is equal to 1, all modes exponentially disappear if 
£iv < 7r/2. This means that the Z 2 symmetric solution (p s = is stable. In 
this case £jv is too small to allow for appearance of the half-compacton + 
half-anticompacton pair at the ends of the chain. 

Let us remark here that such clear cut division into the two possibilities 
is an artifact of the approximation (10). As shown in in the case the 
compacton is obtained from Eq.(8) without this approximation, its width 
depends on O and it approaches tc when (f> — > (the difference vanishes as 
4>q). Therefore, in the real system £at larger than 7r/2 still can be too small 
to accomodate the half-compactons at the ends. In such intermediate case 
certain unstable modes f n ( T ), ^ < n < n , are present but in the final static 
and stable configuration the pendulums are tilted towards one of the rods 
but do not reach it. 

The value of £at can be regulated by changing k. It turns out that evo- 
lution of the system can be investigated analytically also when k is time- 
dependent, k = k(t). This is particularly interesting when k is changed from 
a large value Ki > k c to a small value Kf < k c , where 



is the critical value of the elasticity coefficient at which the first unstable 
mode appears and the symmetric state looses its stability. 

First, we return to the physical position variable x because it is inde- 
pendent of see the first of formulas (6) - equation (28) is written in the 



AmgRaN 2 
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following form 

mgR 

where x G [— Na, Na]. Let us write (f>(x, r) in the form of Fourier series 



d T <f) - -^-d 2 x <j) -0 = 0, (35) 



4>{x,r) = Kij) sin — (x + Na). (36) 
Equation (35) implies that 

h„(r) = \ n h n {r), (37) 

where 

n 2 n 2 K _ n 2 K 

K=1 ~ imgRaW = 1 " ~' (38) 

the dot denotes d/dr. Equations (37) can be integrated also when A n = 
A n (r): 

^n(r) = h n (0) exp (jT A n (<r)da) . (39) 

Let consider in detail the transition which begins at the moment r = 0, 
and such that 

, x I Ki — KT if < T < Tf ,. n . 

«(r) = { - - ; 40 

I K f if r > r/, 

where K > is the constant rate of change of the elasticity coefficient k. 
Hence, 

Ki - Kf 

Tf = : -. 

: K 

In this case formula (39) gives 

K{t) = h n (0) exp (t - —n 2 r + ^-nV ] (41) 



when t < Tf, while at later times 

h n (r) = MO) exp {- {Kl ' Kf l n2 \ exp (r - ^n 2 r) . (42) 
\ 2k c K ) \ k c J 

Formula (41) implies that n-th mode begins to grow at r = r n , where 



However, one has to add the restriction r n < Tf because for larger times 
formula (41) is not relevant. This restriction implies that in fact only modes 
with n < ri/i become unstable, where 

= ^r, < 44 > 

Because Kf > and Ki > k c , 

n/i < 1. 

The term with n — 1 in formula (36) describes pendulums which are tilted 
to one of the rods only - this mode does not give = at any point except 
the ends. Therefore, this mode can not lead to appearance of the kinks. 
Hence, right during the transition, ie. when r < Tf, no compactons and 
anti-compactons are produced. 

During the second stage of evolution of the system, when r > r/, formula 
(42) should be used. It predicts that the number of growing modes remains 
constant, and that it is bounded from above by Uf 2 , where 



n f2 = 



K c 



It is clear that n/2 > 1, and that ri/2 becomes arbitrarily large if Kf is smaller 
and smaller. Below we argue that n/ 2 gives the upper bound on the produced 
compacton - anticompacton pairs. Notice that it does not depend on Ki. Also 
the magnitude of K, or equivalently 77, is not relevant in this respect. 

According to formula (36) each mode evolves independently. Therefore, 
the presence or absence of a given mode is determined by the initial data 
specifying h n (0). If only a single growing mode h n is present, its amplitude 
grows monotonically until the pendulums touch one of the rods. If the friction 
is very large, one may expect that they will not be able to bounce back 
over the energy barrier at = and to fall on the other rod. Therefore, 
in such over-damped case the total number of produced compactons and 
anti-compactons is equal to n — 1. This number does not include the half- 
(anti)compactons which appear at the ends of the chain. The number of 
compactons may differ by 1 from the number of anti-compactons - in such 
the total topological charge is made equal to zero (which is its initial 
value) by the two halves at the ends. For example, if a single compacton is 
present, at both ends there are two half-anticompactons. Thus, in the over- 
damped case the total number of produced compactons and anti-compactons 
is not larger that Uf 2 — I, and it does not depend neither on the initial value 
of k nor on the rate K of the transition. The actual number of produced 
kinks is determined by the initial data h n (0). 
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Formula (41) shows that the value of h n (rf) increases with increasing 
value of the rate K of the transition if n < n c , and decreases if n > n c , where 

2 ^Kf 2 

n r = — n f2 . 

It is clear that n c < nj 2 . Because evolution of the modes at times larger than 
tj does not depend on K, this means that the unstable modes with n < n c 
have larger amplitudes if the transition is rapid. Then, these modes reach 
the angles ±0o earlier than in the case of a slower transition with the same 
initial perturbation. The opposite happens for the modes with n > n c . Thus, 
the magnitude of K has influence on the time of arrival at the final state 
after the transition. 

4 Remarks 

1. One can develop another approach to studying the waves propagating in 
the ground state O (Section 2) by first smoothing the potential V((p) at <fi = 
0o, then applying the standard approach, and removing the regularization in 
the end. The essential difficulty with this approach is that in order to obtain 
waves of non-vanishing amplitude after returning to the original potential 
V{4>) one has to consider waves with finite amplitudes in the regularized 
potential, and such waves probe nonlinearity of the regularized potential. 
Therefore, this approach is in fact more cumbersome than the one adopted 
in Section 2. 

2. It is tempting to compare our results on production of the kinks with 
the well-known phenomenological predictions for the number of produced 
topological defects during phase transitions in condensed matter systems, see 
[TJ [12] for reviews and references, where a strong dependence on the rate of 
the transition is expected, and actually seen in various numerical simulations 
and actual experiments. On the other hand, our explicit calculations of 
that number give results which do not show such dependence. This may be 
explained by the fact that we have tackled only the case of very special initial 
perturbations: the ones with a single eigen-mode present. Clearly, studies 
of this aspect of the model should be continued by considering more generic 
initial states. 

3. Another topic which deserves a separate study is interaction of com- 
pactons, e.g., their scattering. Preliminary investigations reveal a rich dy- 
namics which includes emission of radiation and breather-like states. One 
more interesting aspect not discussed here is interaction of the waves and 
the kinks with the boundaries at ±£at. Recently a study of such interactions 
in another model has appeared P3|. 



13 



5 Acknowledgements 

I would like to thank Danielle Steer for informing me about the link with 
the brane models, Pawel W§grzyn and Zdobyslaw Swierczyriski for useful 
remarks. 

References 

[1] Y. M. Bunkov and H. Godfrin (Eds.), Topological Defects and the Non- 
Equilibrium Dynamics of Symmetry Breaking Phase Transitions. Kluwer 
Academic Publishers, 2000. 

[2] T. W. B. Kibble and M. B. Hindmarsh, Rep. Prog. Phys. 58, 447 (1995). 

[3] H. Arodz, Acta Phys. Pol. B33, 1241 (2002). (Preprint |nTi"n-ps/020100l| ) 

[4] P. Rosenau and J. M. Hyman, Phys. Rev. Lett. 70, 564 (1993). 

[5] F. Cooper, H. Shepard and P. Sodano, Phys. Rev. E48, 4027 (1993). 

[6] A. Khare and F. Cooper, Phys. Rev. E48, 4843 (1993). 

[7] D. Steer, private communication. 

[8] R. K. Dodd, J. C. Eilbeck, J. D. Gibbon and H. C. Morris Solitons and 
Nonlinear Wave Equations. Academic Press, Inc., 1984. Chap. 7. 

[9] S. Dusuel, P. Michaux and M. Remoissenet, Phys. Rev. E57, 2320 (1998). 

[10] T.W. B. Kibble, J. Phys. A9, 1387 (1976). 

[11] W. H. Zurek, Nature, 317, 505 (1985); 

[12] W. H. Zurek, Phys. Rep. 276, 177 (1996). 

[13] N. D. Antunes, E. J. Copeland, M. Hindmarsh and A. Lukas, preprint 
|hep-th/0310103[ 



14 



